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Running of the U{1) conpling in the dark sector 
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The “dark photon” 7^ of a gauged U{l)d can become practically invisible if it primarily decays into 
light states from a dark sector. We point out that, in such scenarios, the running of the U{l)d “fine 
structure constant” otd, with momentum transfer if', can be significant and potentially measurable. 
The 7d kinetic mixing parameter is also expected to run with q^, through its dependence on ad- 
We show how the combined running of e^ad may provide a probe of the spectrum of dark particles 
and, for ad f few x 0.1, substantially modify predictions for “beam dump” or other intense source 
experiments. These features are demonstrated in simple models that contain light dark matter and 
a scalar that breaks U(l)d- We also discuss theoretic considerations, regarding the U(l)d model in 
the ultraviolet regime, that may suggest the infrared upper bound ad ^ 0.1. 


In recent years, various phenomenological considera¬ 
tions have motivated the introduction of a light vector 
boson, the “dark photon” 7 ^, with a mass m-y^ ^ 1 GeV 
- associated with a spontaneously broken gauged U(l)d- 
The dark photon couples to the Standard Model (SM) 
only through kinetic mixing Q of dark charge with hy¬ 
percharge, parameterized by e <C 1. At low energies, 7 ^ 
kinetically mixes with the photon and couples primarily 
to the SM electromagnetic current with a strength ee, 
where e is the electromagnetic coupling. A great deal of 
theoretical and experimental effort has been directed to¬ 
wards dark photon physics |^, Important motivation 
for these efforts have been provided by potential astro- 
physical signals of dark matter Q], as well as the 3.Go- 
deviation from the SM prediction of the muon anoma¬ 
lous magnetic dipole moment 5 ^ — 2 B which can be 
explained by a light jd and e~2xl0“^^. Such small 
values of s most naturally arise from loop effects ||^ of 
particles charged under both the dark and SM hyper¬ 
charge U(l) interactions. For example, a typical I-loop 
value might be £ egd/(167r^). 

If jd is the lightest state in the dark sector then its 
decay is mainly to charged SM states - typically leptons 
£ = e, fi - which could be used to detect it. However, 
in the presence of light dark states, in particular dark 
matter, jd would decay mainly into those particles, since 
in typical scenarios gd ^ £e, where gd is the coupling 
constant of the U{l)d gauge interactions. The different 
phenomenology of such a nearly “invisible” ■jd provides 
new possibilities to explain various anomalies, such as the 
aforementioned g^ —2, but is subject to different sets of 
experimental constraints . In addition, the invisible 
7 d scenario offers an interesting opportunity for produc¬ 
ing and detecting sub-GeV dark matter in accelerator 
based experiments. Boosted dark photons - produced in 
high intensity fixed target experiments - decay in flight 
and lead to a “dark matter beam” which can be detected 


downstream § I |ll|-|l^. The detection rate depends 
on 7 d couplings to the dark and visible sectors, ga and 
ee, respectively; see also Ref. [^. The signal will look 
like electron scattering, modulo radiative corrections and 
scaled by fg^jf. 

The coupling of light dark matter to the dark pho¬ 
ton can provide a mechanism for obtaining its observed 
abundance as a thermal relic density. Here, the annihi¬ 
lation of dark matter via 7 ^ into SM states is mediated 
by kinetic mixing which allows one to derive a rough re¬ 
lation among the model parameters. Using the results of 
Refs. one finds 
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where ad = g^/( 47 r), rrid is the mass of a dark matter 
state, w 10 for a complex scalar ||l^, and w ^ 1 for a 
fermion Here, and elsewhere in this work, ad and 
£ are defined by their values at low momentum transfer 
Given that some of the existing bounds |T^ 
and proposed experiments probe values of e as low as 
10-^ Eq. (0) would then require ^ 1 , keeping other 
parameters at their above reference values, to get the 
correct dark matter relic abundance. Note that if dark 
matter density is set by an asymmetry, efficient annihi¬ 
lation of its symmetric population would require some¬ 
what larger annihilation cross sections, compared to the 
^ pb implied by Eq. (|^), and hence even larger ad val¬ 
ues. Thus, generally speaking, values of aa 0.01 — 0.1 
can be motivated if the dark photon is assumed to decay 
primarily into dark matter states. 

In this work, we point out that the running of ad - 
as a function of momentum transfer q, due to quantum 
loops of light dark sector states - can have important 
phenomenological implications.^ In dark photon models 
it is generally assumed that the kinetic mixing parame¬ 
ter £ corresponds to 0 and remains constant with 
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^ See, e.g., Refs. [ ^ , ^ for possible cosmological effects of run¬ 
ning dark sector couplings, within different frameworks. 
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FIG. 1: Vacuum polarization correction to kinetic mixing. 
The UiVjd current couples to the SM electromagnetic cur¬ 
rent Jgjn, via kinetic mixing, mediated by quantum effects of 
heavy particles Fd,Y (thick loop) that carry both dark charge 
and hypercharge. The vacuum polarization correction from 
light dark states Xd (thin loop) leads to the running of ad{q) 
with momentum transfer q which induces a similar running 
in the kinetic mixing parameter £^(q) 


increasing q. However, as mentioned earlier, kinetic mix¬ 
ing is naturally loop-induced and hence oc aaa, with 
a = e^/(47r). Thus, the running of ad, due to the ef¬ 
fect of light dark states, induces a similar running in 
that would lead to its growth with increasing q and the 
combination ad£^ actually grows like a^- 

The above effect is illustrated in Fig.^, where the vac¬ 
uum polarization correction to the interaction of a U{l)d 
current with the SM electromagnetic current , me¬ 
diated by kinetic mixing, is illustrated. In Fig.^, the thick 
loop, denoted by Fd.v, represents heavy states - charged 
under both the SM hypercharge and U{l)d - that gener¬ 
ate kinetic mixing at the quantum level. The thin loop, 
denoted by Xd, represents light dark states whose vac¬ 
uum polarization contributions lead to the running of 
ad{q)- The loop-induced e^(q) oc a{q)ad{q) inherits its 
running from ad{q)- The well-known running of a{q) 
due to vacuum polarization effects in quantum electro¬ 
dynamics is relatively small and, hence, ignored in our 
discussions. 

Given the discussion following Eq. (||) , values of ad not 
far below unity are well motivated if light dark matter 
communicates with the SM through a heavier invisible 
7d mediator. We show that the growth of aa with mo¬ 
mentum transfer q can constrain the regime of validity 
of calculations, depending on the choice of model param¬ 
eters, and may potentially lead to significant observable 
effects. 

Although the existence of a Landau pole signals strong 
coupling behavior, the onset of strong coupling, where 
perturbation theory starts to break down, is probably 
closer to ttd ~ 1 . From a theoretical point of view, the ap¬ 
proach of a Landau pole for ad can be interpreted as the 
onset of a new non-Abelian interaction, or some other ul¬ 
traviolet completion, that would supplant the low energy 
U{l)d and avoid an ever-growing interaction strength 
with larger q. Such considerations can have interesting 
implications for the underlying model at energy scales 
well above 0, 0 , but will not be studied here, 
except to mention their use as a potential constraint on 
Otd- 


Before going further, we would like to remark that light 
vector models with direct coupling to the SM have also 
been invoked in various contexts Such direct cou¬ 

plings have to be tiny, or else one would observe severe de¬ 
viations from standard physics. Much of the phenomeno¬ 
logical discussion in our work would also apply to these 
models, as long as one assumes 0(1) values for Ud and 
dark sector charges ~ 1 , but small charges for SM 
fields Q™ < 10 “^, under the new Abelian gauge interac¬ 
tions (M . We will assume the alternative kinetic mixing 
picture in our work, as it naturally yields suppressed in¬ 
teractions between the visible and the dark sectors (in 
principle, 7^ could also have mass mixing with the SM 
Z, leading to additional phenomenology [p 0 |). 

As our basic model, we will assume that the dark sec¬ 
tor contains a dark matter state, a fermion "0 or a scalar 
0 , as well as a dark Higgs particle with non-zero vacuum 
expectation value that is responsible for the break¬ 
ing of 17 ( 1 )£;; for simplicity all these particles are assumed 
to have unit charges \Qd\ = 1 under U{l)d- We will fo¬ 
cus on the regime of momentum transfer q > , where 

symmetry breaking effects are negligible and the running 
of ad is significant. In typical proposed fixed target ex¬ 
periments, 10 MeV ^ q ^ GeV, where > 10 MeV 
is probed. As we are focusing on “invisible” dark pho¬ 
ton models, the dark matter state (0 or 0) will be as¬ 
sumed lighter than 7^. Assuming that is not strongly 
self coupled, it is quite natural to expect that its mass 
parameter 771$^ ~ and hence typically less than q. 
Hence, we will include both the dark matter state and 
contributions to the running of ad until q ^ The 

infrared value of the U{l)d coupling can thus be defined 
by its value at q^ = denoted by ad{mj^). 

Since we are focused on the kinematic regime where 
q > rrij^, we will ignore the mass of the vector boson, 
capturing the leading behavior as a function of m^j^/q. 
This suffices for the purposes of our discussions and to 
highlight the key features of generic invisible 7^ scenarios. 
Detailed calculations for specific experiments and model 
parameters lie outside the scope of this work. 

In determining the regime where ^ 1 , higher or¬ 
der effects can become important and we will therefore 
perform a 2 -loop analysis. The 2 -loop beta function of 
U{l)d, with np fermions and ns scalars of unit charge, 
is given by (see, for example, Refs. BH) 








(np + ns) 


( 2 ) 


where Piad) = ^idad/d^i. Here, /x is the renormalization 
scale, which we will later take to be set by the momentum 
transfer q characterizing the interactions of 7^. 

In Fig.|, we have plotted the running of ad in the 
basic model with one dark matter state. Throughout 
our analysis, a dark Higgs scalar assumed to break 
U{l)d, is included in the running and hence 715 > 1 for 
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FIG. 2: Running of ad(q) in the basic model with one dark 
matter state and one dark Higgs boson. The solid (dashed) 
curves correspond to a fermion (scalar) dark matter state and 
the thin (thick) curves correspond to ad{qo) ~ 0.6 (0.9), where 
qo = 0.1 GeV. We have implicitly assumed m.y^ < go- 



FIG. 3: Running of ad{q) with two dark matter states and 
one dark Higgs boson. The solid (dashed) curve corresponds 
to fermionic (scalar) dark matter and ad{qo) = 0.4, where 
qo ~ 0.1 GeV, again assuming < qo. 

all of our results. We will consider a minimum momen¬ 
tum transfer qq = 100 MeV, which is used to set the 
value of ad at the lower kinematic range of typical ex¬ 
periments. The solid and dashed curves correspond to 
fermionic and scalar dark matter [(ni?,ns) = (1,1) and 
(0,2) in Eq. (||)], respectively; ad{qo) = 0.6 (0.9) is rep¬ 
resented by the thin (thick) curve. Based on the form 
of Eq. (^) , for ad^TT the perturbative analysis becomes 
unreliable, due to strong coupling. 

As can be seen from the plot in Fig.H, in all cases, ex¬ 
cept for the one with ad{qo) = 0.6 and {np, ns) = (0, 2), 
the value of ad grows large, i.e. ad ^ tt, by the time q 
reaches ^ a few GeV. This is the typical kinematic do¬ 
main of the proposed fixed target or beam dump exper¬ 
iments. These results suggest that values of ad ^ 0.6, 
used to illustrate the phenomenology in some studies 
MM can typically lead to unreliable predictions, un¬ 
less those values correspond only to the upper kinematic 


range q ^ 1 GeV. In general, numerical predictions are 
more stable when only dark scalars are present in the low 
energy theory. However, even in that case, corresponding 
to (np,ns) = (0,2), we see that the change in ad is not 
negligible for ad{qo) ^ 0.6. 

In Fig.|[ we also present the results for running of ad 
in an extended model that has two dark matter states. 
See, for example. Ref. @ where such models have been 
discussed as a viable setup for sub-GeV dark matter. As 
can be seen from the figure, perturbative analysis be¬ 
comes unreliable for ad{qo) ^ 0.4, with go = 100 MeV, if 
dark matter states are fermionic. However, 2 scalar dark 
matter states do not lead to a significant loss of pertur¬ 
bative validity. Again, the running over q < few GeV 
is not negligible, ^ 50%, even for the all-scalar case and 
can potentially have measurable effects. Overall, we see 
that the assumption of a constant ad over the kinematic 
range of typical experiments can lead to underestimation 
of the predicted rates, even when ad is not close to unity. 

The preceding results point to an interesting possibility 
at fixed target or beam dump experiments, where a range 
of values for momentum transfer q are accessible. For 
Cid{qo) ^ 0.2 —0.3, measurements of processes at different 
q can probe the running of ad{q), in combination with 
the induced running of e{q). As illustrated above, the 
running of modest-sized ad is sensitive to the number 
and type - i.e. fermion versus scalar - of dark states, as 
long as they are well below the typical scale of momentum 
transfer in the measured processes. 

Dehnitive statements regarding the measurement of 
the running with in the above scenarios depend on 
the specifics, such as the production and detection pro¬ 
cesses, and the energy spectrum of the light dark mat¬ 
ter produced by a beam dump or another intense source. 
Nonetheless, key general features of the dark matter scat¬ 
tering through 7 d exchange can be used to outline a po¬ 
tential path toward such measurements, as explained be¬ 
low. Of course, much more detailed studies based on 
specific theoretical and experimental parameters would 
be warranted in designing a dedicated experiment, or 
upon discovery of light dark matter models considered 
here. 

Let us consider the on-shell production of whose 
rate is proportional to e^(m-y^), a constant set by g^ = 
Here, jd primarily decays into dark matter which 
then scatters in the detector with a cross section (Tdm oc 
ad{q)s^{q), where g has a distribution over some range. 
In typical scenarios, kinetic mixing is loop-induced and 
e^(g) oc ad{q), which gives ctdm oc a^(g). 

At g > niry^ , dark matter interactions with the nucleus 
are akin to electromagnetic interactions of a charged lep¬ 
ton with the nucleus, governed by quantum electrody¬ 
namics (QED). This correspondence can be used to ob¬ 
tain precise predictions for the dark matter scattering 
cross section udm- In particular, one could normalize 
(Tdm to the electron (or muon) electromagnetic cross sec- 
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FIG. 4: The running of R/^ as a function of momentum 
transfer q, assuming (A) one and (B) two light fermionic dark 
matter states, corresponding to the solid and dashed curves, 
respectively. We have set ad{qo) ~ 0.25, go = 0.1 GeV, and 
~ 9o is assumed. A scalar (dark Higgs boson) is included 
in the running for both cases. 

tion (Tem oc 1 /g^, which is theoretically well-understood 
and can be precisely measured. Thus, modulo QED ra¬ 
diative corrections and non-zero propagator effects, 
we have (for q > m-yj) 

.R = o’dm/o'em — Od , (3) 

where ^ is approximately a constant (we have ignored 
the negligible running of the QED coupling a over the 
range of considered here). With the above assump¬ 
tions, the value of ^ can be known with good accuracy, 
given the input parameters e and aa at a reference value 
of momentum transfer, which we can naturally take to be 
q = m-y^. In turn, for a given value of one could predict 
the scattering cross section udm and its dependence on 
q^ quite well, in our assumed framework. 

The combined effect of OLd{q) and e^{q) running oc 
Q;^(g) can be significant. For an illustrative numerical 
example, let us consider the case with adiqo) = 0.25 
(assuming ^ go) and go = 100 MeV. With these 
parameters, the effect of running with g^ on the event 
rate can be significant, depending on the value of the f3 
function in Eq. (H). To show this, we consider the cases 
of (A) one dark fermion and (B) two dark fermions con¬ 
tributing to the running (a dark Higgs scalar is included 
in each case). We illustrate the running of i?/^ = a^, for 
the two cases A and B above in Fig.^. One can see that 
the running is significant, over the range g G [0.1,4] GeV, 
in both cases: R/^ changes by a factor of 2 for case 
A (solid curve) and a factor of ^ 4 for case B (dashed 
curve). Furthermore, if sufficient statistics are available 
one can easily distinguish between the two cases, as the 
figure shows, which can potentially probe the dark sector 
spectrum over the measured g^. 

Here, we would like to add a comment. The rise of the 


ratio R with g only encodes the relative increase in (Tdm 
compared to the case with constant couplings. However, 
(Tdm falls like 1 /g^, for g > rrij^ and modulo ad running, 
and hence it is expected that the dark matter scattering 
signal would be stronger for lower values of g^, whereas 
potential backgrounds from neutrino-nucleus scattering 
become more suppressed. The optimal range of g^ for 
detecting the dark matter scattering signals depends on 
the details of the experimental setup. However, as long as 
that range is moderately broad, the g^ running could in 
principle be measurable, as implied by the plot in Figj^. 

So far, we have limited our discussion to the running of 
ad over the GeV-scale values of g, relevant to predictions 
for proposed fixed target (beam dump) experiments. In 
practice, one may not worry if ad becomes too large and 
approaches a Landau pole at g* 1 GeV, as far as those 
predictions are concerned. However, we will argue be¬ 
low that non-perturbative values of ad should typically 
be postponed to g* above the weak scale. To see this, 
note that a straightforward way to resolve the problem 
of an ever-growing ad is to assume the appearance of a 
new non-Abelian gauge interaction Gd at g > g*, whose 
breaking yields U{l)d at lower energies. In typical sce¬ 
narios, the kinetic mixing parameter e vanishes at g*, as 
required by gauge invariance However, to observe 

any events in fixed target experiments, often e ^ 10 “^ is 
required. This implies that e must run to non-zero val¬ 
ues below g*, due to the effects of new states, denoted by 
Fd^Y in Fig.|^, which generate £ 7 ^ 0 at the quantum level 
0. These new states must be charged under both U{l)d 
and hypercharge U{1)y, which means that they cannot 
be lighter than ~ 100 GeV, or else they would have been 
discovered in high energy experiments . 

In light of the above theoretical consistency conditions 
for £ 7 ^ 0 at low energies, we find it well-motivated to 
require that ad should remain perturbative up to g* ^ 
100 GeV. This requirement implies an upper bound on 
the low energy value of the U{l)d coupling ad{qo)- One 
can derive an estimate of this upper bound, using the 
1 -loop running equation 

/ X ^ _ adjq*) _ ... 

1 + ■^ad{q*){nF + ns/4) ln(g*/go)' 

The value of ad{qo) becomes insensitive to the high scale 
value ad{q*) ^ 1 , the onset of strong coupling, as long as 
ln(g*/go) > 1 , and we get 

(2nF-I-ns/2) ln(g*/go) 

In the above, we have implicitly assumed that m.y^ < go 
and we will consider, as before, that go = 0.1 GeV, a 
typical value in the lower kinematic range for fixed target 
experiments. 

We see that for g* = 100 GeV, as a minimal require¬ 
ment based on our preceding discussion, an upper bound 
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otdiqo) ^ 0.68/(ni7’ + 715 / 4 ) is obtained. Hence, for mod¬ 
erate values of 715 and ns one can obtain interesting up¬ 
per bounds. For instance, if 715 = 1 and 715 = 1 (cor¬ 
responding to a dark Higgs), we find ad{qo) ^ 0.5; this 
upper bound gets reduced to ad(qo) ^ 0.3 for np = 2. 
One may entertain much larger values of q*, potentially 
near the Planck scale Mp si 1.2 x 10^® GeV. This could 
be motivated if one assumes that there is no new physical 
mass scale above the electroweak scale, which may pos¬ 
sibly address the stability of the SM Higgs mass against 
large quantum corrections |^, ^ . In that case, one gets 
o^diqo) ^ 0.1/(71f+?^s/ 4) which could place stringent up¬ 
per bounds on adiqo), in the scenarios considered in this 
work. In particular, a dark matter interpretation of the 
'Yd invisible final state implies the relation among model 
parameters given in Eq. (P. For £ < lO”'^, accessible to 
future experiments, that relation would typically require 
Cid ^ 0.1, which would be in perturbative tension with 
q* Mp. 

In passing, we observe that based on constraints from 
Cosmic Microwave Background Radiation (CMBR) mea¬ 
surements 1^, p-wave annihilation of GeV-scale thermal 
relic dark matter is favored, since it becomes less efficient 
at smaller velocities characteristic of the CMBR decou¬ 
pling era 0 . This requirement points to C/(l)ci-charged 
light scalar states as good thermal relic dark matter can¬ 
didates 10, 0]. However, excess annihilation that could 
distort CMBR could also be avoided if dark matter den¬ 
sity is given by an asymmetry at late times, which would 
preclude the possibility of conjugate-pair annihilation. 
We would also like to add that in the models considered 
in our work, dark matter scattering from itself or other 
dark states, mediated by 7 ^, becomes stronger at earlier 
epochs with higher temperatures and larger characteris¬ 
tic values of q. In general, the running of dark sector 
couplings could lead to potentially interesting effects in 
early universe cosmology 00 , but we will not further 
speculate on this question here. 

To summarize, we have considered the running of a 
dark sector U{l)d fine structure constant ad due to the 
presence of light dark particles in quantum loop correc¬ 
tions. Invisible dark photons, assumed to decay on-shell 
to dark matter states, belong to this class of models. We 
observed that the running of ad{q) can lead to a run¬ 
ning kinetic mixing parameter £^( 9 ), for it is naturally 
loop-induced and therefore proportional to ad{q). Un¬ 
der some minimal well-motivated assumptions about the 
dark sector content, we find that, roughly speaking, val¬ 
ues of ttd ^ 0.4 can lead to departure from a perturbative 
analysis and unreliable predictions over the range of mo¬ 
mentum transfer 10 MeV;< 9^1 GeV in proposed fixed 
target or beam dump experiments. We note that ad val¬ 
ues not much smaller than 1 are typical in scenarios 
with sub-GeV dark matter particles that are lighter than 
the dark photon. Light fermionic dark states lead to 
faster running of ad compared to scalar states and can 


result in loss of perturbative control over larger regions 
of parameter space. 

We pointed out that the dependence of ad on q can 
be used to help probe the low-lying spectrum of the dark 
sector in fixed target or beam dump experiments. This 
can be done if measurement of the event rates as a func¬ 
tion of q is feasible, over a moderately broad range of q^. 
Those rates probe the combined running of ad{q) and 
£^( 9 ) oc ad{q) with momentum transfer, leading to a po¬ 
tentially significant sensitivity oc a^{q). We showed that 
for infrared values ad{qo) ^ 0 . 2 , one could expect signif¬ 
icant effects on the event rate from the running, in typi¬ 
cal scenarios. A thorough discussion of q^ running mea¬ 
surements would require input from specific experimen¬ 
tal parameters and would certainly be warranted upon 
the discovery of a light dark matter signal from an in¬ 
tense source. However, on general grounds, we discussed 
how electromagnetic cross sections for electron (or muon) 
scattering from the nucleus can be used to obtain precise 
predictions for the corresponding dark matter scattering 
cross section, in the kinetic mixing scenario considered. 

We also argued that theoretic considerations imply the 
perturbative range for ad should extend to values of q at 
or above the weak scale, perhaps even the Planck scale, in 
which case tighter upper bounds on the low energy value 
of ttd and the dark sector spectrum can be obtained. 

Work supported by the US Department of Energy un¬ 
der Grant Contract DE-SC0012704. 
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